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ABSTRACT 

A successful satellite mission is predicated upon the proper maintenance of the 
spacecraft’s orbit and attitude. One requirement for planning and predicting the 
orbit and attitude is the accurate estimation of the propellant remaining onboard the 
spacecraft. For geosynchronous satellites, a precise propellant-remaining estimation 
is of particular importance. Twenty kilograms (kg) of propellant can add a year to the 
operational lifetime of a satellite such as the Geostationary Operational Environment 
Satellite (GOES)-I. Moreover, the geosynchronous ring is becoming cluttered wit 
propellant-depleted satellites; therefore, an extra 3 to 4 kg of fuel may be required to 
deorbit an expiring satellite out of the geosynchronous ring. For GOES-I, which is 
loaded with over 670 kg of oxidizer and over 420 kg of fuel, accounting for 20 kg of 
propellant requires accuracy in propellant-remaining estimation of within 2 percent. 
Budgeting for the 3 kg of propellant at the end of the mission requires a method with 
an accuracy of within 0*5 percent. 

This paper focuses on the three methods that were developed for calculating the pro- 
pellant budget: in particular, the errors associated with each method and the uncer- 
tainties in the variables required to determine the propellant remaining that 
contribute to these errors. Based on these findings, a strategy will be developed for 
improved propellant-remaining estimation. The first method is based on Boyle’s law, 
which relates the values of pressure, volume, and temperature (PVT) of an ideal gas. 
The PVT method is used for both the monopropellant and the bipropellant engines. 
The second method is based on the engine performance tests, which provide data that 
relate thrust and specific impulse (I sp ) associated with a propellant tank to that 
tank’s pressure. Two curves representing thrust and specific impulse as functions o 
pressure are then generated using a polynomial fit on the engine performance data. 
The third method involves a computer simulation of the propellant system. The pro- 
pellant flow is modeled by creating a conceptual model of the propulsion system con- 
figuration, taking into account such factors as the propellant and pressurant tan 
characteristics, thruster functionality, and piping layout. 

Finally, this paper presents a thrust calibration technique that uses differential cor- 
rection with the computer simulation method of propellant-remaining modeling. 
Thrust calibration will provide a better assessment of thruster performance and 
therefore enable a more accurate estimation of propellant consumed during a given 
maneuver. 
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1. METHODS OF COMPUTING PROPELLANT REMAINING 


Introduction 


Before presenting the detailed descriptions of the propellant estimation methods, a brief introduction 
on the propulsion system operation during a spacecraft mission is in order. 

For the liquid propellant/gas pressurant system considered in this paper two modes of operations are 
feasible: blowdown and pressure regulated. In blowdown mode, the propellant tank is pressurized by 
the gas pressurant and then during the mission the propellant tank pressure is allowed to decay as 
propellant is "blown out” of the tank when the thrusters are firing. In pressure regulated mode, the 
propellant tank pressure is maintained constant by supplying additional pressurant gas into the pro- 
pellant tank during the thruster firing. The advantage of the pressure regulated mode is the constant 
propellant flow rate which is necessary to maintain in a bipropellant type engine to ensure constant 
mixture ratio for optimum thruster performance (see Section I, Method 1). 

The propellant estimation methods described here can be used to model an engine operating in both 
blowdown and pressure regulated mode, blowdown mode only, and pressure regulated mode only 
(Methods 1, 2, and 3 respectively). 

Two phases of the mission are mentioned in this paper: the transfer orbit phase— when the satellite is 
maneuvered to achieve mission orbit; and the station keeping phase— when the satellite is maneuvered 
to maintain the mission orbit. For GOES-I, 86 percent of the propellant is used during the transfer 
orbit phase (NASA phase). GRO on the other hand is inserted into the mission orbit by the launch 
vehicle; therefore, ideally, 100 percent of propellant is used for station keeping and controlled reentry. 

Method 1: Pressure, Volume, Temperature 

The PVT method is based on an assumption of ideal gas behavior of the pressurant gas. Since the 
pressurant gas is helium, the ideal gas approximation is valid. Boyle’s law is then used to estimate the 
propellant remaining based on the amount of pressurant that was forced into the propellant tank. The 
procedure is as follows: 

Suppose that the volume of propellant displaced from the tank is equal to the volume of pressurant 

forced into the tank, given that pressurant and propellant do not mix and that the tank volume does not 
change. Then, 



dM f = rhof • dVHe.f 

(1-1) 

where 

dMf = propellant forced out of the tank 

rhof = propellant density 

dVfie.f = volume of pressurant forced into the tank 


Using Boyle’s law to compute the pressurant volume and writing the propellant density as a function of 
tank temperature, we have 


dy u [ _ R • dM He • Z He • T He 

^He 

(1-2) 

where 

R = pressurant (helium) gas constant 

dm He = mass of pressurant forced into the propellant tank 

^He - pressurant compressibility 

Tpie - pressurant temperature 

PHe = pressurant pressure 
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id 


rhof — bo + bi'T + b2'T 2 


(1-3) 


iere t = propellant tank temperature 

bo, bi, b2 = propellant density coefficients 

nee pressurant and propellant are in the same tank and do not mix. 


id 


T He = T 


P He = P ' ? 


sat 


(1-4) 

(1-5) 


here P = propellant tank pressure 

p sat = propellant saturation pressure 

ropellant saturation pressure may be expressed as a function of temperature, T, 


(1-6) 


Psat = 10 (a °" ai/T 32//T2) 

here ao, ai, a 2 = propellant saturation pressure coefficients 

ressurant compressibility may be expressed as a function of pressurant pressure, P He , and tempera- 

rn 


ire, T He , 


Z H e = 1 + P H e ’ [P + r ' (T - T s )] 


(1-7) 


'here /J, r = pressurant compressibility coefficients 
T s = pressurant standard temperature 

"he expression for propellant used may then be rewritten as a function of tank pressure and tempera- 
jre and the mass of pressurant forced into the propellant tank, 


dM f = R • dM He • [b 0 + bi ■ T + b 2 • T 2 ] • 
+ ^- T + r T 2 -r*T s -T 


p _ IQ ^ 30 ” ai /^ ~ ^P 1 ) 


(1-8) 


Hie accuracy of the propellant estimation using this method is only as good as the certainty in the 
values of P, if and dM He - The uncertainty in the tank pressure and temperature is based on the ratings 
jfthTpressure transducers and the temperature sensors as well as the 

uncertainty in the pressurant mass forced into the propellant tank can be consi ® , 

as the uncertainty in the pressurant mass leaving the pressurant tank (assuming t ere we n 
which is a function of the Lading conditions and the telemetry readings of the pressure and tempera- 
ture of the pressurant tank. 


Since pressurant gas behaves as an ideal gas, Boyle’s law applies as follows: 


V 0 


ITlr 


ZHe 


PHe ’ VHe 
mHe ’ Z He ’ T He 


(1-9) 
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where P 0 

T 0 

V 0 

m 0 

^He,o 

P He 

T He 

V H e 

m He 

^He 


= loading pressurant tank pressure 
= loading pressurant tank temperature 
= loading pressurant tank volume 
= , loading pressurant mass 
= loading pressurant compressibility 
= pressurant tank pressure 
= pressurant- tank temperature 
= pressurant tank volume 
= pressurant mass 
= pressurant compressibility 


Using the fact that pressurant tank volume is a function of pressurant tank pressure and temperature, 


V 0 = V h + qi • (P 0 - P h ) + q 2 • (T 0 - T h ) (1-10) 

V He = V h + qi • (P He - P h ) + q 2 ■ (T He - T h ) ( 1 - 11 ) 


where 


Vh = standard pressurant tank volume 
Ph — standard pressurant tank pressure 
Th = standard pressurant tank temperature 
qi = pressure coefficient 
q 2 = temperature coefficient 


and the fact that compressibility is also a function of tank pressure and temperature, Equation (1-7), 
the amount of pressurant forced into the propellant tank may be expressed as follows: 


dM He = m 0 - m He 


( 1 - 12 ) 


or 


dM He = mo • 


1 


p He • Tp Vh + qi • ( p He ~ P h) + q2 ’ (T He - T h ) 
p 0 • T He Vh + qi • (Po - Ph) + q 2 • (T 0 - T h ) 


1 + P 0 • \p + r • (To - T s )] 

1 + P He • \p + r • (T He - T s )) 


(1-13) 


The error in propellant estimation may be expressed in terms of the uncertainties in propellant tank 
temperatures and pressures obtained from telemetry and the uncertainty in the pressurant mass 
forced into the propellant tank. Therefore, the error in propellant estimation due to these variables 
can be defined in standard fashion as 


_ <5 (dMf) . 

5 (P) Uncertainty 

(l-14a) 

_ <5 (dM f ) , 

6 (T) * uncertainty 

(1- 14b) 
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<5 (dMf) ,w H 

Error dM He = ^ (dM He ) dMH uncertainty 


(l- 14 c) 


The uncertainty in pressurant mass can be expressed as la error in the pressurant mass due to the 
uncertainties in the pressurant tank loading conditions and telemetry uncertainty in pressurant t 
pressure and temperature: 


dMH ep ncertainty 


= /(ErrorpHe ) 2 + (Error THe ) 2 + (Error Po ) 2 + (Error To ) 2 + (Error mo ) J (l- 14 d) 


where 


Error-pHe - 


Errorp 0 = 


Errorjo = 


Error mo = 


• P He uncerta inty 

<5 (P H e) 

(l-14e) 

^^-4 • THe unc ertainty 

(l-14f) 

o (P 0 ) 

(l-14g) 

6 (dMne) . 

1 Ouncertainty 

0 (T 0 ) 

(l-14h) 

d (dMne) . _ 

mouncertainty 

(1-MO 


6 (mo) 


Blowdown mode operation may be simulated by assuming that there is no change in pressurant mass 
in the propellant tank (i.e., pressurant tank is shut off). Then, the propellant remaining becomes a 
function of the propellant tank pressure and temperature change where dM He - constant. 

In a bipropellant propulsion system such as the one in GOES-I, the pressurant forced into a given tank 
is a function of the split ratio (the ratio of pressurant mass forced into the two tanks) as weU as 
function of the pressurant mass leaving the pressurant tank. The split ratio is, m turn, a function o e 

pressures and temperatures of the two propellant tanks and the mixture ratio ( the ratio ot ^ the mass 
flow rates of the two propellants). The pressures and the temperatures are obtained from telemetry 
the mixture ratio is defined by the manufacturer to ensure the optimum thruster performance. The 
expression for the split ratio in terms of the above quantities is derived as follows. 

„ ,. „ . dMnei 
Split Ratio = 


(1-15) 


dMne2 

where dMnei = mass of pressurant forced into propellant tank 1 
d M H e2 = mass pressurant forced into propellant tank 2 

Combining Equations (1-1) and (1-2) and introducing subscripts to distinguish between the two pro- 
pellant tanks, 

ZHel ' T Hel 


= rh0fl • * 

dM He i 


PHel 


(l-16a) 
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and 


dM c 

dM He 2 


rh 0f2 • — ' Zlie2 ’ THe2 
PHe2 


(l-16b) 


Then, using the definition of the mixture ratio (MR) and Equations (1-3) and (1-7) for propellant densi- 
ties and pressurant compressibility, respectively, 


Split Ratio = MR • _^gl + ’ T2 + bp • T2 2 _ PHei . T He 2 

b^ + b 0 i • T1 + b o2 * Tl 2 P He2 T He i 

(1-17) 

. I ± PHe2 • H + T ■ (T2 - T s )l 

1 + PHel • W + r • (Tl - T s )] 

where pressurant partial pressure can be expressed in terms of the pressure and temperature of the 
propellant tank using Equations (1-5) and (1-6). Then the pressurant forced into either propellant tank 
as a function of the total pressurant leaving the pressurant tank and the split ratio is 


dM He i 


dM He 


Split Ratio 
(Split Ratio + 1) 


(l-18a) 


dM He 2 = dM He 


1 

(Split Ratio + 1) 


(1- 18b) 


Therefore, the set of error equations (Equation (1-14)) for a given propellant tank (1 or 2) must be 
expanded to include the errors in pressurant mass forced into the propellant tank due to the uncertain- 
ties in the mixture ratio and the uncertainties in 
in the system. That is, 



Errorp! 

_ <5(dM He ) . P1 

r ^uncertainty 

(l-19a) 

Errorxi 

_ <5(dM Hc ) . 

A A uncerta j nt y 

(l-19b) 

Errorp2 

_ <5(dM„ e ) . 

Sf 2 \ ^uncertainty 

(1-1 9c) 

Errors 

_ d(dM H e) . _ 

^ ^uncertainty 

(1-1 9d) 

Error MR 

(5(dM He ) 

<5(MR) MR uncertain ty 

(l-19e) 
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lie 1(7 error in pressurant forced into a given propellant tank then becomes 

dMHe^ruin* = [(Error P „ e ) 2 + (Error ™) 2 + (Error Po ) 2 + (Error To ) 2 + (Error ^) 2 ^ 

+ (Error P i ) 2 + (Error T i ) 2 + (Errorp 2 ) 2 + (Error T2 ) 2 + (Error MR ) 2 ] 1/2 

Jsing the data for the GOES-I propellant system operating in pressure regulated mode as ^ example: 
iven nominal operating propellant tank pressures and temperatures (indexes 1 and 2 indicate 
xidizer and the fuel tanks, respectively) of 

PI = P2 = 230 psi T1 = T2 = 20° C 

nd pressurant tank loading mass, pressure, and temperature 


P 0 = 3300 psi 


T 0 = 21° C 


m c = 2.54 lbm 


nd assuming the operating pressurant tank temperature remains constant 

T He = T 0 

he pressurant tank pressure at the end of NASA phase is reduced to 

P He = 200 psi 

nd the optimum mixture ratio as supplied by the manufacturer is 

MR = 1.610 

\ssuming that there are no uncertainties in loading conditions of the pressurant tank and mass, the 
lartials are computed to be 


= -1.4 lbm/psi 

<5(P2) 


= 0.57 lbm/° K 

5(T2) 


<5(dMf) _ non C 

<5(dMHe2) 


The partials to compute the uncertainty in pressurant mass forced into the fuel tank are 


<S(dMHe2 ) _ _q. 0004 lbm/psi 

<5(P H e) 

<3(dM He 2) = _ 0 002 lbm/psi 

m) 

<3(dM H e2) _ +0.002 lbm/psi 

<5(P 2 ) 

^ dMHa) = -0.3 lbm 
d(MR) 


<3(dM He 2) = + 0 003 lbm /° k 

6 ( T He ) 

^( dMHe2 I = +0.002 lbm/ 0 K 
<5(Ti) 

= -0.001 lbm/° K 

<5(T 2 ) 
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1 1 ,li 


The uncertainties in pressure and temperature due to the telemetiy resolution are 

Pluncertainty(telem) = ± 0.4 psi Tl uncertainty(te , em) = ± 1.0° K 

P2 U ncertainty(telem) = ± 0.4 psi T2 uncertainty(teIetn) = ± 1,0“ K 

PHe unC ertainty(telem) = ± 0.4 psi THe U ncertainty(telem) = ± 1.0° K 

The uncertainty in pressure due to the transducer accuracy and the resulting total root-mean-square 
uncertainty in pressure readings are M 


Pluncertainty(trans) ““ ± 

1.35 psi 

Pluncertainty(tot) — 

± 1,41 psi 

- 

P2 uncer ^ a i n jy^j rans ^ = i 

1.35 psi 

P2 uncerta] - nt y( tot ) = 

± 1.41 psi 

- 

P He uncerta j nty ( trans ) = ± 

1.35 psi 

PHe uncerta j nty ( tot ) = 

± 1.41 psi 

j 


The uncertainty in the mixture ratio as supplied by the manufacturer is 

MRyncertajpLty = ± 0.024 


Thus, the errors in fuel used associated with the resulting uncertainties in the fuel tank pressure tem- 
perature, and amount of pressurant forced into the fuel tank ( 3o error in pressurant forced into the 
fuel tank is ±0.027 Ibm) are * 

Errorp 2 = 2.0 Ibm Error T2 = 0.57 Ibm Error dM He = 9.17 lbm = 

This shows that the fuel-used estimate for the GOES-I spacecraft in the pressure-regulated mode has a . 

3a uncertainty of ± 28.2 lbm. Thus, given a 911.6 lbm estimated fuel usage, the relative error in fuel- ■ 
remaining estimation is 3.1 percent. 

Method 2: Thrust and Specific Impulse Performance Data 

This section presents the mathematical argument for the thrust and specific impulse (I SD ) curves 
method of computing the propellant consumed from a tank during a specified time interval. These 5 

curves describe thrust and I sp as functions of pressure and temperature. The method assumes that the 
propellant system behaves according to the thrust and I sp performance curves derived through empir- E 
ical testing of the propellant system. The equations to describe these curves are derived through poly- 
nomial fitting and take on the following form when the first three terms of the polynomial are used: = 

+ c4 • P 1 

(1-20) S 


d3 + d4 • P 

(1-21) * 


F = cO + cl • P - c2 • P 2 • 


Tref 


cl 


I S p = d0 + dl • P - d2 • P 2 


Tref 
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vhere F = thrust 

: specific impulse 

= propellant tank pressure 
= inlet propellant temperature 

= inlet propellant temperature at which the data were taken 
- thrust polynomial coefficients 
= I sp polynomial coefficients 

Note that if T and Tref are equal, then thrust and I sp are functions of tank pressure only.) 

jiven thrust and I sp , the propellant flow rate is easily determined: 


Isp 

P 

T 

Tref 

cO, cl, c2, c3, c4 
dO, dl, d2, d3, d4 


Q) = F/L 


( 1 - 22 ) 


vhere d> = propellant flow rate. 

rhe propellant mass escaping the tank during a certain time period is 


ti 

dMf = 

to 


a) • dt 


(1-23) 


where dt = time period. 

Substituting Equations (1-20) and (1-21) into Equation (1-22), integrating with respect to pressure and 
temperature, dividing by the change in pressure and temperature to obtain the average flow rate, and 
then substituting into Equation (1-23), we find that 


dM f = 


ti 


to 


r T1 

PI wdT 

Jl° dP 

TI - TO 


f P1 

. P0 


PI - P0 


dt 


(1-24) 


where P t o = tank pressure at tO 
P tl = tank pressure at tl 
T t o = tank temperature at tO 
T t i = tank temperature at tl 

Note that when there is no change in temperature or pressure in a given time interval, the flow rate is 
COTStant with respect to that variable over this time interval, and therefore, the integration step with 
respect to the unchanging variable should be omitted. 

Since thrust and I SD method depends on tank pressure variation, it is meaningful to use this method 
on"“n .he blowdoSn mode of operation when a significant change in tank pressures can be obsetved. 

The errors associated with this method are inherent to the instruments used in d e rivingthrastandl,p 
data points as well as the data regularity required to produce a close polynomial fit. Assuming that 
thrust and l sp are well-behaved functions and that the instruments used to take the data are extreme y 
accurate, the error in determining the propellant flow rate is then a function of un certainties n t e 
burn start and stop time and of the uncertainties in tank pressure and temperature at burn start and 

stop time. 
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Errorpo 


Errorp! 


Errorxo 


Errorji 


Error t0 


Error t i 


- <5 (dMf) . 

6 (P0) ru “ n “rtamty 

(1-25 a) 

_ a ( dM f) . p. 

ft /pj\ ^uncertainty 

(l-25b) 

_ 15 ( dM f) . Tn 

6 (TO) 1 ^uncertainty 

(l-26a) 

_ 6 (dMf) . 

5 (Tl) ' ^uncertainty 

(l-26b) 

_ d (dMf) . , n 

6 (tO) ^uncertainty 

(l-27a) 

_ d (dMf) . ,, 

^ uncertainty 

(l-27b) 


The following example for the Gamma Ray Observatory (GRO) satellite, which operates in blowdown 
mode, shows the calculations of the fuel used from a tank given that tank’s pressures at start and end of 
burn and the burn duration. 

Using GRO main satellite thrusters performance coefficients and assuming that the tank temperature 
remains the same as the reference temperature during which the curves data were taken: 

F = (3.3502 Ibf) + (0.39898 lbf/psi) • P - (0.0001463 Ibf/psi 2 ) • P 2 


I sp = (222.52 s • g) + (0.064329 s • g) • P - (0.0000672 s ■ g) • P 2 

where the units of thrust and I S p coefficients are as appropriate, and, given a 2-minute ascent maneuver 
and propellant tank pressures at the start and end of the maneuver of 

dt = 120 sec 
P0 = 400 psi 
PI = 334 psi 

the partials are computed to be 


(5(dMf) 

= 0.072 Ibm/psi 

(5(dM f ) 

= 0.548 lbm/sec 

<5(P0) 

m 

d(dMf) 

a(Pi) 

= 0.071 lbm/psi 

<5(dMf) 

= -0.548 lbm/sec 
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The uncertainties in pressure and time due to the telemetry resolution are 

POuncertainty " ±3.0 psi tO U ncertainty = ± 0.256 Sec 

Pluncertainty = ±3.0 psi tl U ncertainty = ± 0.256 sec 

This shows that the fuel-used estimate from a GRO tank during an ascent maneuver has a 3cr uncer- 
tainty of ± 1.1 lbm. Thus, computing that the total fuel used during the ascent from that tank is 
65 8 lbm the relative error in the fuel-remaining estimation is 1.7 percent. Note that this method is an 
approximation that relies on engine performance to follow the curves obtained during ground testing. 

Method 3: Conceptual Model of the Propulsion System 

A conceptual model of the propulsion system involves creating a schematic representing the layout of 
the propellant piping, tank, and thruster configuration. Then a set of mathematical expressions must 
be developed to describe the physics of this system, using the data obtained from the manufacturers on 
such system characteristics as the flow resistance through the piping, characteristic propellant velocity 
and thrust coefficients for all thrusters, the throat areas of the thrusters, and the temperature and 
pressure of the tanks. A good example of the development of such a model is the GOES-I bipropellant 

system model. 

The GOES-I propellant system consists of a pressurant tank, a fuel tank, an oxidizer tank, one main 
satellite thruster (MST), and 12 attitude and orbit control thrusters (AOCT) arranged in strings A and 
B, each containing six AOCTs. 

A model representation of the bipropellant system consisting of only one thruster (e.g., the MST) may 
be used to derive the following set of governing equations representing the physics of the system (Fig- 
ure 1): 


Pc = Pf - Z K f"f 

(l-28a) 

Pc = Po _ 

(l-28b) 

F = A t CfP c 

(l-28c) 

F = (a ) 0 + (Of) Isp 

(l-28d) 

F = A + B • I S p 

(l-28e) 


where V Ko = line resistances (oxidizer) 

V Kf = line resistances (fuel) 

At = throat area 

A, B = coefficients in Equation (l-28e) 
P Q = oxidizer tank pressure 

Cf = orifice coefficient 

P f = fuel tank pressure 
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Figure 1 . Schematic of the Simplified Bipropellant System Containing Only 
the Main Satellite Thruster 


d) { 

= fuel flow rate 

E 


■ oxidizer flow rate 

= 

Pc 

= chamber pressure 

1 

F 

= thrust 

* 

!sp 

= specific impulse 

= 


This simple model can then be expanded to include the entire system. The conceptual model repre- 
senting the system functionalities is shown in Figure 2. = 

Using the conceptual model in conjunction with the propellant system’s physical constants, a set of 
governing equations relating propellant flow rates and thruster chamber pressure can be derived for 
each of the 13 thrusters in the same manner as for MST r 

Then Ihe mathematical representation of the functionality of the whole system is accomplished in 
combining the above equations for AOCTs and MST by applying a physical constraint of propellant 
flow continuity inherent to the system. That is, propellant mass flowing into a junction is equal to 

propellant mass flowing out of that junction. For example 

<bft = 6){ + 0)fA + <i>fB (1*29) 

Solving the system of equations described above will give the propellant flow rates through each ^ 

thruster and the chamber pressure of each thruster. Propellant used due to each thruster is then the t 

product of the flow rate and thruster on time. The propellant remaining may also be calculated by 
using the chamber pressure in the thrust and I S p performance data curves for each thruster. 
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Figure 2. Conceptual Schematic of the GOES-I Propulsion System 


Currently, this model is used for propellant estimation in pressure regulated mode of operation. How- 
ever, by solving the system of equations for each new pressure reading in the propellant tanks, this 
model may be used for propellant estimation in blowdown mode. 

The uncertainty in the conceptual model method comes mainly from the error in the flow resistance 
and thruster coefficients, tank pressure and temperature, and thruster on time. Also, there is added 
error in any method chosen to solve the system of nonlinear equations. 

The error due to the method of solving the equations is simply the smallest tolerances of the variable 
under which the method converges to a solution. The error due to flow resistance and thruster 
coefficients, tank pressure, and temperature is determined by adding maximum error to these vari- 
ables and then solving the equations to see the amount by which the solution under maximum error 
deviates from the nominal solution obtained by using nominal values of these variables. The error due 
to time uncertainty is simply the product of the time uncertainty and the computed propellant flow 

rate. 

The following is an example of error in GOES-I MST firing in pressure regulated mode propellant- 
used prediction as computed by the bipropellant engine model using nominal propellant flow 
resistance and thrust coefficients as supplied by SS/Loral and assuming nominal tank pressures and 
temperatures 

Ko = 110.713 lbf • s 2 /lbm • in 5 K f = 190.944 lbf • s 2 /lbm • in 5 


C F = 1-865 P 0 = Pf = 230 psi T 0 = T f = 21.3° C 
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and based on the following uncertainties: uncertainty in fuel and oxidizer tank pressures due to trans- 
ducer accuracy and telemetry resolution (total pressure uncertainty is root-mean-squared of these 
two); uncertainty in temperatures due to telemetry resolution, 

P uncertainty(telem) = ± 0.4 psi T unC ertainty = ± 1.0* K 

Puncertainty(trans) = ± 1.35 psi P U ncertainty(tot) = ± 1.41 psi 

and uncertainty in MST propellant flow resistance (Ko and K f ) and thrust (C F ) coefficients as given by 
SS/Loral (Reference 1) 

K°uncertainty = ± 0.541 Ibf * S 2 /lbm • jn 5 

^uncertainty = ± 2.583 Ibf • s 2 /lbm • in 5 
CFuncertainty = i 0.00236 

The resulting fuel and oxidizer flow rates and la errors in the flow rates due to the above uncertainties 

combined with 0.00001 convergence tolerance of the flow rates when solved for using Runge-Kutta 
method are 


< = 0.21857 Ibm/s Error^ = 0.00090 lbm/s 

<i)f = 0.13566 lbm/s Error^ = 0.00087 lbm/s 

The uncertainty in MST on time due to the telemetry resolution is 

^uncertainty = i 0.023 Sec 

Assummgthenominal GOES-I first two apogee maneuver firings, the total MST on time is 96 minutes 
( 5,760 sec). Then the fuel and oxidizer masses used (as computed by the bipropellant engine model) 

dM 0 = 0.21857 lbm/s • 5760 s = 1258.96 Ibm 

dM f = 0.13566 lbm/s - 5760 s = 781.40 lbm 

The errors in the propellant-used computations are: 
from errors in flow rates: 

Erro W„ - ± 5-18 tt>»> 

Error dMf,», = ± 5 01 lbm 

from errors in thruster on time: 

Error dMo ,t = 0.005 lbm 

Error d Mf,t = 0.003 lbm 


6130-29 


494 



Thus, based on the total propellant used during the burn, the 3 o error in fuel used is ± 15.03 lbm, the 
3a error in oxidizer used is ± 15.54 lbm, the relative error in the fuel-used estimation is 1.9 percent, 
and the relative error in the oxidizer-used estimate is 1.2 percent. 


2. CALIBRATION TECHNIQUES 


The models discussed in Part 1 of this document neglect to take advantage of the actual performance 
data of the spacecraft during the mission, which can be determined from the actual orbit achieved 
after the burn or the status of the orbit during the burn. That is, actual thrust delivered by the engines 
can be deduced from the orbital data available through tracking. This section presents a method of 
calibrating the conceptual propellant system model by using the actual thrust of the GOES-I satellite 
as determined from the orbit achieved. 

Biprop differential corrector is a FORTRAN program that modifies the parameters of the GOES 
Bipropellant engine model (developed using the algorithm described in Part 1; see Reference 2) until 
the solution for total thrust obtained by the model matches the observed total thrust produced by the 
GOES propulsion system. The correction applied to the parameters is based on the information, 
according to SS/Loral, that the engine components most likely to vary during a burn are the propellant 
flow resistance coefficients for the MST section of the piping. Since the propellant piping is such that 
there are no isolated thrusters, the AOCTs are also affected by the varying MST resistance coefficients. 
However, the AOCTs are not affected when the MST is off, because MST off indicates zero propellant 
flow to the piping with varying resistance coefficients. Thus, Biprop differential corrector is used only 
when the MST is on — that is, during the NASA phase of the mission. In summary, the Biprop differen- 
tial corrector is designed to correct for total thrust produced by the MST and the AOCTs combina- 
tions by adjusting the propellant (both oxidizer and fuel) flow resistance coefficients of the MST 
piping. The single constraint on varying the MST fuel and oxidizer resistance coefficients, given by 
SS/Loral, is that the mixture ratio (the ratio of the oxidizer flow rate and the fuel flow rate) for the MST 
must equal a predetermined constant. This section discusses (1) calculation of total thrust and average 
MST mixture ratio taking into account the AOCT duty cycles; (2) differential corrector requirements, 
(3) the differential corrector algorithm; and (4) some examples to illustrate the function and perform- 
ance of the Biprop differential corrector program. 

2.1 CALCULATION OF TOTAL AVERAGE THRUST AND AVERAGE MST 
MIXTURE RATIO 

During a burn, the AOCTs are usually fired for a shorter time period than the MST. The on time of the 
AOCTs is described by a duty cycle (percentage of the burn time that the AOCTs are on). The MST 
stays on for the entire burn period. The equation for the total thrust is the average of all thrusters that 
are firing weighted according to each thruster’s on time. As was shown in the study of the effects of 
multiple thruster firing on thruster performance (see Reference 3), for the total thrust magnitude cal- 
culations it is valid to assume an average duty cycle for all the AOCTs that are on. Likewise, it is valid to 
assume that all AOCTs that are on are firing at the same time and at the beginning of the burn. There- 
fore, the equation for total weighted average thrust is a sum of two parts: one for the MST firing alone 
and another for the MST firing together with the AOCTs. Hence, 


T av = T M ST(oR) ' (1 “ Dut y Cycle/100) + (T M sT(on) + X T aoc) 


• (Duty Cycle/100) 


(2-1) 
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where T av 

T MST(off) 

T MST(on) 

X Taoc 

Duty Cycle 


= total weighted average thrust 
= MST thrust while AOCTs are off 
= MST thrust while AOCTs are on 

= total thrust produced by AOCTs 

= percent of the burn time that the AOCTs are on 


The mixture ratio used in the differential corrector must also be averaged, taking into account the duty 
cycles. The mixture ratio in the MST changes when the AOCTs go on, because the flow rates to the 
MST are changed. Therefore, the average MST mixture ratio must be calculated in a fashion similar to 
average, total thrust calculations. The formula for the weighted average MST mixture ratio is 


MRav — MR.MST(off) (1 - Duty Cycle/ 100 ) + MRmst(oii) ’ (Duty Cycle/ 100 ) ( 2 - 2 ) 

where MR av = average weighted MST mixture ratio 

MRMST(off) = MST mixture ratio while AOCTs are off 
MRMST(on) = MST mixture ratio while AOCTS are on 

This average mixture ratio is constrained to equal a predetermined value. 

2.2 DIFFERENTIAL CORRECTOR REQUIREMENTS 

The differential corrector algorithm has two requirements for the function that describes the system: 

1. The function must be continuous over a chosen interval. 

2. The function must be differentiable on this interval. 

Both of these requirements must be true for total thrust and MST mixture ratio as functions of the 
MST flow resistance coefficients. Since the system being modeled is a physical system, the thrust pro- 
duced by the system must be directly related to the propellant flow in the system. From the governing 
equations of the bipropellant engine model (see Reference 4) we have for any given thruster 


- p f - X (K© ‘ ®f(i) 2 ) 

(2-3) 

— P© - ^ (Ko(i) " ^o(i) 2 ) 

(2-4) 

T = (A, • C f ) • P c 

(2-5) 


where P c 

= 

chamber pressure 

Pf 

= 

fuel pressure 

Po 

= 

oxidizer pressure 

K f 

= 

fuel resistance coefficient 

Ko 

= 

oxidizer resistance coefficient 

d) { 

= 

fuel flow rate 



oxidizer flow rate 

T 

= 

thrust 

A t 

= 

throat area of the thruster 

c f 

ss 

thrust coefficient 
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’herefore, from Equations (2-3) and (2-4), P c (K f ) and P c (Ko) are linear functions. Since T is propor- 
ional to P c , T(Kf) and T(K<j) are also linear functions. Thus, taking Kf(MST>and Ko(MST) to be variable 
esistances for the MST, the functions that relate these resistances to the total thrust, T(K f (MST)) and 
’(Ko(mst))> are linear by the above argument, and thereby meet the requirements of the differential 
orrector algorithm. 

ntuitively, there must be a smooth relationship between the flow resistance coefficients and the flow 
ates. That is, the flow rate of a propellant in a pipe is smoothly related to the resistance of the pipe’s 
nterior surface. The mixture ratio is simply the oxidizer flow rate divided by the fuel flow rate. Also, 
he flow rates are never zero, since the MST is always on when the differential corrector is required, 
fherefore, the mixture ratio is a smooth function of the oxidizer and fuel resistance coefficients, and 
hus meets the requirements of the differential corrector algorithm. 

1.3 DIFFERENTIAL CORRECTOR ALGORITHM 

[Tie problem of thrust correction is defined by two variables and two constraints. The two variables are 
he fuel resistance coefficient and the oxidizer resistance coefficient. The two constraints are that the 
:omputed thrust must equal the actual thrust and that the computed mixture ratio must equal the 
ictual mixture ratio. The requirement that mixture ratio be fixed implies that the the actual mixture 
atio is equal to the nominal mixture ratio within a specified tolerance. The following procedure must 
>e used in performing the differential correction on the bipropellant engine model: 

1 Obtain the actual mixture ratio (MRa) using bipropellant engine model with nominal flow 
resistance coefficients supplied by the manufacturer. 

2. Obtain the actual, total thrust magnitude from the calibrated maneuver mode. 

3. Until the computed total thrust (T n ) and the computed mixture ratio (MR n ) are within the 
specified tolerance of the constraints, iterate with 

Kf(MST) (i + 1) = Kf(MSTXi) + AKf(MST) (2-6) 

Ko(MST)(i+l) = Ko(MSTXi) + AK<j(MST) (2-7) 

,yhere AKf^yr) and AK^mst) are obtained via the differential corrector method 

_-l r 

AKq(mst) _ 6T/ dKo^MST) <5T/(5Kf(MST) T a - T n 

AKf(MST) <5MR/<5 Ko(mst) <5MR/ 6Kf(MST) MR a - MR n 

where <5 K<j(mst> <5Kf(MST) = perturbation applied to the coefficients 

<5T, <5MR = (Tpert - Tn), (MRpert - MR n ) respectively 

Tpert, MRpert = thrust and mixture ratio, respectively^ computed by BIPROP 

using perturbed coefficients Kf(MSTXi) + 6Kf(MST). 

Ko(MSTXi) + 6Ko(MST) 

T n , MR n = thrust and mixture ratio, respectively, computed by 

BIPROP using unperturbed coefficients 
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2,4 EXAMPLES OF BIPROP DIFFERENTIAL CORRECTOR 

Table 1 gives some examples of the differential corrector performance. The numbers used for ti 
actual thrust (T a ) were chosen only for the testing purposes. The nominal resistance coefficients we 
part of the data given by the manufacturer for testing the bipropellant engine model software. / 
shown in the table, the differential correction makes good progress in only two to three iterations wit 
relative error for thrust and mixture ratio specified at under 0.5 percent. Note that in Case 2 we fin 
significantly higher resistance than in Case 1, although the difference between actual and total thru, 
in both cases is almost the same. This may be understood as follows; 

Let 


DTmst - T M ST(off) - T M ST(on) (2-9) 

so the DTmst is the change in the MST thrust caused by AOCTs firing. Then, substituting DT^.t int 
Equation (2-1), we get 5 

Tav = T MST (off) + (Duty Cycle/100) ■ (£t A oc - DT MS t) (2-10) 


Table 1 . Differential Corrector Performance 


Thrusters 

Firing 

Nominal 
resistanc 
coeff fcie 
as suppli 
FACC 

e 

nts 
ed by 

Nominal 
total 
thrust 
computed 
by BIPROP 

Actual 

average 

mixture 

ratio 

computed 

by 8 1 PROP 

Actual 

total 

thrust 

from 

cal i bra ted 

manuever 

mode 

Actual 
resistance 
coef f icients 
computed by 
Differential 
Corrector to 
obtain actual 
thrust 

Total 
number of 
i t era t ions 
performed by 
Differntial 
Corrector to 
obtain the 
actual thrusi 

names 

Oxidizer 
^o(MST) ^ 

Fuel 

<lc f <MST) } 

< T n> 


<v 

Oxidizer 

^o(HST)^ 

Fuel 

<IC f(MST) } 

n 

MST only 

110,713 

190.944 

103.876 

1.611 

93.6 

176.079 

297.618 

3 

MST, 2A 

110.713 

190.944 

110.748 

1.644 

95.0 

178.624 

303.305 

2 

MST, west 
face AOCs 

110.713 

190.944 

112.613 

1.643 

102,0 

141.416 

255.606 

2 

MST, all 
AOCs 

110.713 

190.944 

116.034 

1.642 

140.0 

53.046 

95.709 

3 


NOTE ' For the examples here, the differential corrector uses the following inputs; 

Specified Relative Error for T a and MR a : 0.5% 

Duty Cycle: 10% 

Bum Time; 3,000 sec 
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Although DTmst > 0, we observe that Taoc - DTmst) > 0 • Also, from Equations (2-3), (2-4), 
and (2-5), we know that T vtsi'(off) decreases with increasing resistances. Likewise, from Equa- 
tion (2-10), we have 

TMST(off) = T av - (Duty Cycle/ 100) • (X ^aoc - DTmst) (2-11) 

This shows TMST(off) < T av in Case 2 , while TMST(off) = T av ' n Case 1. Therefore, we expect Case 2 to 
require higher resistances. Moreover, as the resistances increase, the flow rates to the MST decrease, 
and, from continuity conditions, the flow rates to AOCTs must increase. Thus, 

(Duty Cycle/100) ■ Taoc — DTmst) increases as resistances increase. Similarly, 

(Duty Cycle/100) • (^T T A oc - DT MS t) increases as the number of the AOCTs firing or the Duty 

Cycle increases. That is, the increase in required resistances to obtain the same average thrust be- 
comes more marked as the number, or the duty cycle, of the AOCTs firing with the MST increases. 

Case 4 shows that the resistance coefficients may also be decreased in order to account for a better 
thruster performance than expected. 

When such a calibration technique is used, the propellant flow rates and the chamber pressure of each 
thruster that was firing during the burn are adjusted in the process to reflect more closely the actual 
performance of the thrusters. Using the information of case 1 from Table 1 as an example, if the actual 
thruster performance is 86 percent of the predicted performance than the difference in fuel used com- 
puted using nominal and calibrated flow resistance coefficients is 107 lbm. For GOES-1 107 Ibm is 
equivalent to 2.5 years of mission lifetime. Therefore, the propellant-remaining calculations derived 
by the conceptual model of the system are more realistic, since the model reflects the actual perform- 
ance of the propellant system as observed during the mission. 

REFERENCES 

1. Ford Aerospace, Space Systems Division, SSD-TR 108 17 A, GOES IJK/LM Propulsion Require- 
ments Analysis, DRL 303-01, June 1988 

2. Computer Sciences Corporation, SEAS Quick Note No. GOES-468-89020, Expanding the 
Bipropellant Model to Include 12 AOC Thrusters, prepared by S. Torgovitsky, May 12, 1989 

3. — , SEAS Quick Note No. GOES-468-89006, Effects of Multiple Thruster Firing on Thruster Per- 
formance and Fuel Remaining, prepared by S. Torgovitsky, December 7, 1988 

4. — , CSC/TR-89/6003, Geostationary Operational Environmental Satellites (GOES)-I/M Flight 
Dynamics Support System (FDSS) Volume II.' Functional Specifications, Revision 1, C. Audain 
et al., February 22, 1989 


*U 5 GOVERNMENT PRINTING OFFICE: 


1991. 52?. oe^sooa 


6130-29 


499 


r NASA 

National Aeronautics and Report Docu mentation Paqe 

Space Administration 1 

1 . Report No. 

NASA CP-3123 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Flight Mechanics/Estimation Theory Symposium - 1991 

5. Report Date 
October 1991 

6. Performing Organization Code 
554 

7. Authors ) 

Thomas Stengle, Editor 

8. Performing Organization Report No. 

91B00133 r 

10. Work Unit No. 

9. Performing Organization Name and Address 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 

1 1 . Contract or Grant No. 

13. Type of Report and Period Covered 
Conference Publication 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Adm 
Washington, D.C. 20546-0001 


inistration 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Thomas Stengle is Head, Attitude Analysis Section, Flight Dynamics Analysis Branch, Goddard Space Flight Center 
Greenbelt, Maryland ' 


16. Abstract 

This conference publication includes 26 papers and abstracts presented at the Flight Mechanics/Estimation 
Theory Symposium on May 21-23, 1991. Sponsored by the Flight Dynamics Division of Goddard Space Flight 
Center, this symposium features technical papers on a wide range of issues related to orbit-attitude prediction, 
determination, and control; attitude sensor calibration; attitude dynamics; and orbit decay and maneuver strategy. 
Government, industry, and the academic community participated in the preparation and presentation of these 
papers. 


I 17. Key Words (Suggested by Author(s)) 

18. Distribution Statement 

Flight Mechanics 

Spacecraft Dynamics 

Unclassified 

Estimation Theory 

Orbit Determination 


Attitude Determination 

Mission Analysis 

Subject Category 13 


19. Security Classif. (of this report) 

1 

20. Security Classif. (of this page) 

21. No. of pages 

22. 

Unclassified 

Unclassified 

508 

i 


A22 


NASA FORM 1626 OCT 86 



